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ABSTRACT 

In compact white dwarf (WD) binary systems (with periods ranging from minutes 
to hours), dynamical tides involving the excitation and dissipation of gravity waves 
play a dominant role in determining the physical conditions (such as rotation rate and 
temperature) of the WDs prior to mass transfer or binary merger. We calculate the 
amplitude of the tidally excited gravity waves as a function of the tidal forcing fre- 
quency u = 2{Q — Q S ) (where f2 is the orbital frequency and fl s is the spin frequency) 
for several realistic carbon-oxygen WD models, under the assumption that the outgo- 
ing propagating waves are efficiently dissipated in outer layer of the star by nonlinear 
effects or radiative damping. Unlike main-sequence stars with distinct radiative and 
convection zones, the mechanism of wave excitation in WDs is more complex due to 
the sharp features associated with composition changes inside the WD. In our WD 
models, the gravity waves are launched just below the helium-carbon boundary and 
propagate outwards. We find that the tidal torque on the WD and the related tidal 
energy transfer rate, Etide, depend on oj in an erratic way, with -E t idc varying by orders 
of magnitude over small frequency ranges. On average, £Wc scales approximately as 
f2 5 u; 5 for a large range of tidal frequencies. 

We also study the effects of dynamical tides on the long-term evolution of WD 
binaries prior to mass transfer or merger. Above a critical orbital frequency fi c , corre- 
sponding to an orbital period of order one hours(depending on WD models), dynam- 
ical tides efficiently drive £l s toward 0, although a small, almost constant degree of 
asynchronization (Q — Ct s ~ constant) is maintained even at the smallest binary peri- 
ods. While the orbital decay is always dominated by gravitational radiation, the tidal 
energy transfer can induce significant phase error in the low-frequency gravitational 
waveforms, detectable by the planned LISA project. Tidal dissipation may also lead 
to significant heating of the WD envelope and brightening of the system long before 
binary merger. 
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1 INTRODUCTION 

Compact white dwarf (WD) binary systems (with orbital 
periods in the range of minutes to hours) harbor many inter- 
esting and unanswered astrophysical questions. An increas- 
ing number of such systems are being discovered by recent 
surveys (e.g. Mullally et al. 2009; Kulkarni & van Kerkwijk 
2010; Steinfadt et al. 2010; Kilic et al. 2011; Brown et al. 
2011; see Marsh 2011 for a review). The orbits of these sys- 
tems decay via the emission of gravitational waves, which 
could be detected by the planned Laser Interferometer Space 
Antenna (LISA) (Nelemans 2009). Depending on the WD 
masses and the physics of the merger process, these merging 
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WD systems may produce single helium-rich sdO stars, giant 
stars (R CrB stars), stable mass transfer AM CVn binaries, 
or possibly undcrluminous supernovae. Most importantly, 
compact WD binaries in which the total mass is near the 
Chandrasekhar limit are thought to be the probable progen- 
itors of type la supernovae upon a stellar merger at the end 
of the orbital decay process (Webbink 1984; Iben & Tutukov 
1984). Recent studies have provided support for this "double 
degenerate" scenario (e.g., Gilfanov & Bogdan 2010; Di Ste- 
fano 2010; Maoz et al. 2010) and even sub-Chandrasckhar 
WD mergers may lead to type la supernovae (van Kerkwijk 
et al. 2010). 

Prior to merger, tidal interactions may affect the prop- 
erties of the binary WDs and their evolutions, including the 
phase evolution of the gravitational waves. Previous stud- 
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ies have focused on equilibrium tides (e.g., Iben et al. 1998; 
Willems et al. 2010), corresponding to quasi-static deforma- 
tion of the star. Such equilibrium tides are unlikely to play 
a role in the tidal synchronization/dissipation process. Iben 
et al. (1998) estimated the effect of tidal heating in the WD 
based on the assumption that the (spherically averaged) lo- 
cal heating rate is equal to the rate of rotational energy 
deposition required to maintain synchronization. They sug- 
gested that the binary WDs may brighten by several mag- 
nitudes before merger. 

In fact, in a compact WD binary, as the orbital decay 
rate due to gravitational wave radiation increases rapidly 
with decreasing orbital period, it is not clear if tidal effects 
are sufficiently strong to drive the binary system toward 
synchronous rotation. The critical orbital period for syn- 
chronization is unknown. For this reason, the majority of 
recent WD merger simulations (e.g., Segretain et al. 1997; 
Loren-Aguilar et al. 2009; Pakmor et al. 2010,2011) have 
assumed the merging WDs to be non-synchronized prior to 
merger. However, whether the WDs are spin-synchronized 
may affect the merger product and the possible supernova 
signature: for example, the strong velocity shear between 
the stars upon contact would be significantly reduced for 
the merger of a synchronized binary. The degree of synchro- 
nization also determines the tidal luminosity of the binary 
prior to merger. Indeed, it is possible that tidal dissipation 
contributes significantly to the brightness of some of the re- 
cently observed WD binaries (e.g., Brown et al. 2011). 

In a recent paper (Fuller & Lai 2011, hereafter Paper 
I), we used linear theory to calculate the the excitation of 
discrete gravity modes in a WD due to the tidal gravita- 
tional field of a compact companion star (a WD, neutron 
star or black hole). The existence of discrete modes requires 
that gravity waves be reflected near the surface of the WD. 
In this case, tidal energy and angular momentum transfers 
between the WD and the binary orbit occur only during a se- 
ries of resonances, when the g-mode frequency a a equals 2Q, 
(where Q, is the orbital frequency). Our calculations showed 
that while the dimensionless (mass-weighted) amplitude of 
the resonantly excited g-mode is not extremely non-linear 
(it approaches ~ 0.1), the displacement associated with the 
mode becomes large in the outer layer of the WD where 
the density is low. In other words, while the mode does not 
reach a non-linear amplitude in the bulk interior of the star, 
it becomes very non-linear in the outer layers even before 
resonance. We concluded that tidally excited gravity waves 
are likely to continually damp in the outer layer of the WD, 
preventing the formation of discrete modes. A proper treat- 
ment of dynamical tides in binary WDs must take account 
of this continuous wave damping. 

In this paper, we calculate the tidal excitation of gravity 
waves in binary WDs assuming that the waves are efficiently 
damped in the WD envelope. To this end, we implement an 
outgoing wave boundary condition near the WD surface. 
Unlike gravity modes (which have a set of discrete eigen- 
frequencies), the outgoing wave boundary condition permits 
the excitation of gravity waves at all frequencies, and thus al- 
lows for a continuous process of tidal dissipation. Similar cal- 
culations have been implemented for early-type stars (Zahn 
1975,1977; Goldreich & Nicholson 1989) and solar-type stars 
(Goodman & Dickson 1998; Ogilvie & Lin 2007). In early- 
type stars, gravity waves are excited at the boundary be- 



tween the convective core and radiative envelope, propagate 
outwards and dissipate in the outer envelope. In solar-type 
stars, gravity waves are similarly excited at the interface be- 
tween the radiative core and convective envelope, but prop- 
agate inward before dissipating (via non-linear wave break- 
ing) near the center of the star (Barker & Ogilvie 2010,2011). 
Unlike main-sequence stars, WDs do not contain a simple 
two-zone structure of convective and radiative regions, and 
it is not clear how and to what extent gravity waves are ex- 
cited. The outgoing wave outer boundary condition allows 
us to calculate the rate at which energy and angular mo- 
mentum are transferred to the WD as a function of orbital 
period. We can then calculate the orbital period at which 
tidal effects can compete with orbital decay due to gravita- 
tional radiation. At this orbital period, the synchronization 
process can begin. Furthermore, by scaling our results to 
rotating WDs, we can determine the WD spin period and 
energy dissipation rate at any orbital period. 

This paper is organized as follows. In Sections 2-4 we de- 
rive the basic equations for tidally forced stellar oscillations, 
the boundary conditions and the tidal angular momentum 
and energy transfer rates. In Section 5 we discuss our nu- 
merical method and present several test calculations, where 
we emphasize the importance of using a self-consistent stel- 
lar model in order to obtain reliable amplitudes for tidally 
excited gravity waves. In Section 6 we present our numeri- 
cal calculations of tidal excitations for realistic WD models. 
Since our numerical results reveal a complicated dependence 
of the tidal energy transfer rate on the tidal frequency, we ex- 
amine a simple semi-analytic model in Section 7 to shed light 
on the mechanism of gravity wave excitation. In Section 8 we 
use the results of previous sections to study the long-term 
spin-orbit evolution of WD binaries, including spin synchro- 
nization, the tidal effect on the low-frequency gravitational 
radiation waveforms and tidal heating of the WDs. We con- 
clude in Section 9 with a discussion of theoretical uncertain- 
ties and future work. 



2 BASIC EQUATIONS 

The dynamical tide of the WD (mass M) is driven by the 
external gravitational potential of the the companion (mass 
M'). The leading order (quadrupole) potential is 

U«(r,t) = t/(r)[y 22 (M) e - 1 ^ +i? 2 (MK w *] (1) 

with 

... . GM'W 2 2 2 

U(r) = - z r . (2) 

Here a is the orbital separation, uj — 2Q. is the tidal fre- 
quency for a non-spinning WD (we will account for the spin 
effect in Section [sj) , Q. is the orbital frequency, and W22 = 
\/37r/10. The actual fluid perturbations in the WD can be 
written as £ ac (r, t) = £(r, t) +£* (r, t) for the Lagrangian dis- 
placement and <5P ac (r, i) — SP(r,t) + 5P* (r,t) for the Eu- 
lerian pressure perturbation, and similarly for other quanti- 
ties. In the following, we shall consider perturbations (£, SP, 
etc.) driven by the potential U{v,t) = U(r)Y 22 (8, ^)e" 1 "'. 
We shall adopt the Cowling approximation (so that the grav- 
itational potential perturbation <5<I> associated with the den- 
sity perturbation is neglected, i.e., <5$ = 0), which is valid for 
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gravity waves in the star. We will consider adiabatic oscil- 
lations, for which the Lagrangian perturbations in pressure 
and density are related by AP — a 2 Ap, where a s is the 
adiabatic sound speed. This is a good approximation in the 
bulk of the star where the thermal time is much longer than 
the wave period. 
Letting 



5P(r,t) = 6P(r)Y 22 (6,i 



and 



£(r,i) = [&(r)f + e±(r)rV_ L ]y a3 (M)e- 
the fluid perturbation equations reduce to 



(3) 
(4) 



+ 

a% pa 2 



L <V-fe^ = o, ( 5 ) 



and 



-8P' + -^-5P+(N 2 
p pal 



2 )& + U' = 0, 



(6) 



where the ' denotes d/dr. In equations ([5| and Q, Li and 
N are the Lamb and Brunt- Vaisala frequencies, respectively, 
given by [note we will continue to use the notations Li, 1(1 + 
1), and m, although we focus on I — m — 2 in this paper] 



2 l(l + l)a 2 s 



and 



N 2 = - 2 ' dp 



1 



dP a 2 



(7) 



(8) 



The other perturbation variables are related to 8P and £ r 
by 



5 P =±5P+^ r . 
at 9 

"1/2^2, 



Defining Z — \ r £r, where 
pa 2 s \u) 2 

equations |5| and Q can be combined to yield 
Z" + k 2 (r)Z = V(r). 

Here, 



(9) 
(10) 

(11) 
(12) 



+ 



9 



+ 



and 



V(r)= X - 1/2 



-(g/a 2 sY 
g/al 

Kl + i) f-x 1 



X 9_ 

X "1 



X a 2 J a 2 



U. 



(13) 



(14) 



In the WKB limit \k\ > l/H and fcj > l/r, where 
H — \P/P'\ ~ a 2 /g is the pressure scale height, equation 
( 13 1 simplifies to 



k 2 (r) = 



;{L 2 -LU 2 ){N 2 ~UJ 2 ). 



(15) 



This is the standard WKB dispersion relation for non-radial 
stellar oscillations (e.g., Unno et al. 1989). For d 2 -C L 2 and 



u) 2 <C A^ 2 , the wave equation \12\ reduces to 



z n , l(l + l)(N 2 -u 2 ) z 
r 2 u) 2 



_ 1/2 /(/ + l)iV 2 u 



-x 



(16) 



Then, as long as \Z"/Z\ 3> |x"/xl (which we expect t o be 
true because Z" ~ — k 2 Z and %" ~ x/H 2 ), equation (16 1 



is identical to the oscillation equations used by Zahn (1975) 
and Goodman & Dickson (1998). 



3 BOUNDARY CONDITIONS 



Equations (JsJ and (JsJ) or equation ( 12 I can be solved with 
the appropriate boundary conditions at r — r ou t near the 
stellar surface and at r — n n — > at the center of the star. 



The general solution of equation (121 can be written as 

Z(r) = c+Z+(r) + c-Z-ir) + Z cq (r), (17) 

where c+, c_ are constants. Z+(r) and Z-(r) are two inde- 
pendent solutions of the homogeneous equation Z" + k 2 Z = 
0, and Z Bq (r) represents a particular solution of equation 
(12 1. We choose the outer boundary r out to be in the wave 
zone (k 2 > 0). If k 2 (r) varies slowly such that \k' /k\ <C k > 
0, then the two independent WKB solutions to the homoge- 
neous equation are 



Z±(r) = ~= exp | = / 



kdr 



(18) 



where r is an interior point (r < r out ). For w 2 < i^, the 



WKB wave dispersion relation [equation ( 15 1] reduces to 



uj 2 = N 2 k 2 _/(k 2 + k I L ), where kj_ = /(Z + l)/rV*vhich implies 
that the radial component of the group velocity is —ujk/(k 2 + 
kj_). Thus, withoi > and k > 0, Z- oc e~ l ^ r ° kdr represents 
an outgoing wave, while Z+ oc e > T a kdr represents an ingoing 
wave. An approximate particular solution of equation ( 12 1 



Y__}_(YL 
¥ fe 2 \¥ 



(19) 



where the second term is smaller than the first by a factor of 
ikH) 2 or (kr) 2 . This represents the "non-wave" equilibrium 
solution^ 

Throughout this paper, we adopt the radiative condi- 
tion at the outer boundary (r = r out ), i.e., we require that 
only an outgoing wave exists: 



Z(r) 



Z eq (r) + 



sfk 



exp 



kdr 



(20) 



This implicitly assumes that waves propagating toward the 
WD surface are completely damped by radiative diffusion 
(Zahn 1975) or by non-linear processes. We will check this 



Note that the equilibrium tide usually refers to the f-mode re- 
sponse of the star to the tidal force. Here, we use the term "equi- 
librium" to refer to the "non-wave" solution. 
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assumption a posteriori from our numerical results (see Sec- 
tion 6.5 1. Thus, near the outer boundary, the radial displace- 



ment £ r behaves as (for oj <C L t 



€r(r) 



M 2 



-Z(r) 



■. exp 



^Jpr 2 (N 2 -oj 2 ) 
Here ££ q represents the equilibrium tide 



u 



N' 2 



N 2 



1 - 



2gr 



l(l + l)a 2 N 2 



kdr). (21) 



(22) 



where we have retained only the first term of equation ( 19 1. 
For oj 2 < iV 2 , this further simplifies to £° q ~ -U/g (Zahn 
1975). The constant c_ specifies the amplitude of the out- 
going wave which is eventually dissipated in the stellar en- 
velope; this is the constant we wish to determine from nu- 
merical calculations. 



In practice, to implement equation (211 at the outer 



boundary, we require a very accurate calculation of the non- 
wave solution ££ q . This can become problematic when the 
conditions \k\ 2> 1/H and u <C Li are not well satisfied. 
Since the transverse displacement £x for gravity waves is 
much larger than the radial displacement in the wave zone, 
it is more convenient to use £j_ in our outer boundary con- 
dition. We define 



Mr) 



\ V 2 



(23) 



with D = N 2 — uj 2 . Equations (JHJl and Q can be rearranged 
to yield 



(24) 



where 



kl{r) 



1 

4 

9 



hi 



In - 



pr 
~D 

r 2 
D 



OJ 

+ — + 



In 



pr 



D 

l(l + l)D 



9 



(25) 



and 



Vi(r)=- 



N 



pr 

u 

u 



1/2 



hi 



r 2 N 2 
Dg 



U 



(26) 



For kx > 1/H and oj 2 <C Lf, the functions kj(r) and Vi(r) 
simplify to 



and 



K(r) 



V 1 (r) ~ - 



1(1 + 1){N 2 - oj 2 ) 



pr 



1/2 



D fUr 2 N 



(27) 



(28) 



Again, adopting the radiative boundary condition at r — 
r ou t, we have 



-z , , Vx 1 (Vx 

f\j-\ \ \ 1 



+ ^— exp 
kx 



kxdr), (29) 



where c_ is a constant. Thus, the transverse displacement 
£x(r) behaves as 



D\ 1/2 1 



k± 
pr 2 



Mr) 



1/2 



exp 



kidr 



(30) 



where several constants have been absorbed into c_. The 
equilibrium tidal transverse displacement Cj q ( r ) i s given by 

'Ur 2 iV 2x ' 



l(l + l)r\ g D 

For uj 2 < iV 2 , this reduces to (for Z = 2) 

l/t/r 2 ^' 
6r V 3 



(31) 



(32) 



in agreement with Goldreich & Nicholson (1989). Thus, we 
implement the radiative boundary condition at r — r out as 



pr 2 /kj 



2(pr 2 /fci) 



— ikx 



(33) 



with £x computed from £ r and SP using equation j9j. 

The inner boundary condition can be found by requiring 
the radial displacement to be finite at the center of the star. 
This requires 



I 



\^ + U J (Nearr = 0). 



(34) 



4 ANGULAR MOMENTUM AND ENERGY 
FLUX 

As the wave propagates through the star, it carries an angu- 
lar momentum flux to the outer layers. At any radius within 
the star, the z component of the time-averaged angular mo- 
mentum flux is 



j z (r) = ( j> dQr p(Sv r + 8v*) (8v$ + 5v*^)r sin#y, (35) 

where (...) implies time averaging. With 

8v r = — iio£ r (r)Yi m e~ luJt (36) 



and 



x . t , v _ v -tut mw£ ± (r) _ ibJt 

6v^ = -iw£±(r)rV <j,Yi m e = — — Y lm e 

smO 



(37) 



we find 



j z {r) = 2 I dnr :i puj 2 Re\iC(r)Y; m m^ ± (r)Y lr 



= 2mo; 2 pr 3 Re( 



(38) 

In the wave zone, the fluid displacement consists of 
an equilibrium ("non-wave") component and a dynamical 
(wave) component, £ = £ eq + £ dyn . Since the equilibrium 
tide component is purely real (assuming negligible dissipa- 
tion of the equilibrium tide), Re^i£J: q (,±\ = 0, and the 
equilibrium tide does not contribute to angular momentum 
transfer. The cross terms Refe yn *^ q ) , and Re(^° q *£ ± yn ) 
are opposed by a nearly equal and opposite Reynold's stress 
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term (see Goldreich & Nicholson 1989) and do not con- 
tribute significantly to angular momentum transfer. Thus, 
the Refi^f 3 ™ £^ yn ^ term dominates angular momentum 
transferrl Equation ( 38 1 then becomes 



j z (r) =2ma;V 3 Re[iC, dyn *£x n ]- 



(39) 



In the outer layers of the WD where £ dyn is a pure 
outgoing wave (oc e~ %kr ), equation ^ can be rearranged to 
obtain the relationship between £ dyn and ^f yn in the WKB 
approximation (k 3> 1/H) with lj 2 <C L 2 : 



dyn 



kr 



1(1 + 1) 

Then the angular momentum flux is 



dyn 



(40) 



Jz 



2ml(l + l)^\tr\ 2 

4y/l(l + i)to 3 pr 3 d 2 
TV IS± 1 • 



(41) 



where we have used the dispersion relation (equation 271 
with w 2 <C TV 2 and set m = 2. This expression agrees with 
that found in Goldreich & Nich olso n (f989). From the scal- 
ing of £ dyn provided in equation |33| , we see that the angular 
momentum flux is constant (independent of radius) in the 
outer layers of the star. Since the wave pattern frequency (in 
the inertial frame) is (the orbital frequency), the energy 
flux carried by the wave is given by E = Q.J Z . 

Once we have solved our differential equations |5] and 
[6| with the appropriate boundary conditions, we can use 
equation ( 39 1 to determine where angular momentum and 



energy are added to the wave, i.e., where the wave is gener- 
ated. In the WD interior, the waves travel both inwards and 
outwards and thus carry no net angular momentum, so the 
value of J z oscillates around zero. However, near the outer 
boundary, the value of J z is constant and positive because 
there only exists an outgoing wave. The region where the 
value of J z rises to its constant value is the region of wave 
excitation, because it is in this region where energy and an- 



gular momentum are added to the waves (see Section 6.3| 



The energy and angular momentum carried by the out- 
going wave is deposited in the outer envelope of the star. 
Thus, the constant values of J z and E near the outer bound- 
ary represent the net angular momentum and energy trans- 
fer rates from the orbit to the WD. Since £ dyn oc M' /a 3 , the 
angular momentum and energy transfer rates can be written 
in the form 



where 



Jz = ToF(w), E = T QF(lj) 



(42) 



(43) 



and F is a dimensionless function of the tidal frequency u 
and the internal structure of the star. For WDs with rotation 
rate fi S) the tidal frequency is uj — 2(fi — fi s ). 



2 It can be shown that in the WKB limit the Reynold's stress 
associated with the dynamical response is negligible. 



5 WHITE DWARF MODELS 

Figure [I] depicts three WD models provided by G. Fontaine 
(see Brassard et al. 1991). These WD models are taken from 
an evolutionary sequence of a M = 0.6Mq WD, at effective 
temperatures of T = 8720K, T = 5080K, and T = 2720K. 
The WD has a radius R ~ 8.97 x 10 s cm and a carbon- 
oxygen core surrounded by a 10~ 2 M helium layer, which 
in turn is surrounded by a 10~ 4 M layer of hydrogen. The 
models shown have been slightly altered in order to ensure 



thermodynamic consistency (see Section 6.1 1 



The Brunt- Vaisala frequency can be expressed as 



P Xp 



V + B) 



(44) 



where the symbols have their usual thermodynamic defi- 
nitions, and the Ledoux term B accounts for composition 
gradients (see Brassard et al. 1991). In the core of the WD, 
the value of iV 2 is very small due to the high degeneracy 
pressure, which causes \t in equation ( |44[ | to be small. The 
sharp spikes in TV 2 are created by the carbon-helium and 
helium- hydrogen transitions, and are characteristic features 
of WD models. These sharp features in realistic WDs make 
it difficult to construct toy WD models or to understand 
how gravity waves propagate through the WD. From Figure 
[l] it is evident that cooler WDs have smaller values of N 2 
throughout their interiors. However, the spikes in iV 2 have 
little dependence on WD temperature because they are pro- 
duced by composition gradients rather than thermal gradi- 
ents, thus these features are unlikely to be strongly affected 
by tidal heating. 



6 NUMERICAL CALCULATIONS OF TIDAL 
RESPONSE 

6.1 Numerical Method and Importance of 
Self-Consistent Stellar Model 

To calculate the amplitude of the gravity waves excited in 
a WD by its companion, we integrate the inhomogeneous 
equations Q and Q with the appropriate boundary condi- 
tions given by equations ( 33 1 and ( 34 1. We use the relaxation 
method discussed by Press et al. (2007). The integration re- 
quires a grid of points containing stellar properties (p, iV 2 , 
a 2 , g) as a function of radius, and solves the equations on a 
grid of (possibly identical) relaxation points. 

When creating the grid of data points representing the 
stellar structure, one must be very careful in ensuring that 
the stellar properties are consistent with one another. In 
particular, the Brunt Vaisala frequency is given by 



N = -g 



(45) 



If the value of iV 2 in our stellar grid is not exactly equal 
to the right hand side of the above equation as calculated 
from the values of p, g, and a 2 , the stellar properties will 
not be self-consistent. Such inconsistency may arise from 
the inaccuracy of the stellar grids, or from the interpolation 
of the stellar grids (even if the original grids are exactly self- 
consistent). We have found that even a small inconsistency 
can lead to large error in the computed wave amplitude. The 
reason for this can be understood by examining equation 
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£ 10" 2 

CD 
Q 




0.4 0.6 
Radius (R) 



Figure 1. The square of the Brunt- Vaisala (thin solid lines) and Lamb (dashed lines) frequencies (for / = 2), in units of GM/R? , and 
the density (thick solid line) as a function of normalized radius in three WD models. The models arc taken from an evolutionary sequence 
of a DA WD with M = 0.6M Q , R = 8.97 X 10 3 km, and effective temperatures of T = 8720K (top), T = 5080K (middle), and T = 2720K 
(bottom). The spikes in the Brunt Vaisala frequency are caused by the composition changes from carbon to helium, and from helium to 
hydrogen, respectively. Note the formation of a convective zone just below the carbon-helium transition zone as the WD cools. 



Tidal Synchronization and Dissipation in White Dwarf Binaries 7 



(??). When tracing back to equations (JsJ and d6j, one can 
see that the N 2 term on the right-hand side is actually the 
sum of two terms. That is, the value of N 2 on the right-hand 



side of equation ( 16 1 is calculated via equation (45 1 from our 
grid of p, a 2 , and g values, while the N 2 term on the left 
hand side of the equation is taken directly from our grid of 
iV 2 values. If these two values of iV 2 differ (by the amount 
8N 2 ), then a "false" excitation term will be introduced on 



the right-hand side of equation (16 1, given by 
l(l + l)5N 2 U 



Vf 



(46) 



This false term can vary rapidly with radius depending on 
the error in the stellar grid. In Section [7] we discuss how 
sharp changes in the excitation term can be responsible for 
the excitation of the dynamical component of the tidal re- 
sponse. Thus, the false excitation term introduced by even 
small numerical inconsistencies can cause large errors in cal- 
culations of the dynamical tide. 

To test our methods, we calculated the tidal response 
for a simple massive star model. The results are discussed 
in Appendix [X] and are consistent with previous studies of 
gravity waves in massive stars (e.g., Zahn 1975,1977 and 
Goldreich & Nicholson 1989). 

6.2 Calculation with Toy White Dwarf Model 

To understand wave excitation in WDs, we first examine a 
toy model constructed to mimic the structure of a WD. Ex- 
amining the Teff = 8720K model, we see that it contains a 
sharp rise in iV 2 at the carbon-helium boundary, preceded 
by a small dip in iV 2 near the top of the carbon layer. Conse- 
quently, we have created a toy model with a similar dip and 
rise in TV 2 in the outer part of the star. To create this model, 
we first construct a smooth density profile (identical to that 
of an n — 2 polytrope, along with a smooth N 2 profile that 
mimics the dip-rise features associated with the C-He tran- 
sition in real WDs. Next, we compute a thermodynamically 
consistent sound speed profile using the equation 

_ _ ^ 
\dP/d~p ~ p 5 " 



2 ( 

a, — 



(47) 



Since the density profile is that of a polytrope, the dP/dp 
term can be calculated analytically. 

We solve the forced oscillation equations as a function 
of the tidal frequency uj. Figure [3] shows the energy flux and 
wave amplitude as a function of radius for a given value of 
uj. The small oscillations in energy flux are due to imperfect 
numerical calculation of the dynamical component of the 
wave and do not actually contribute to energy or angular 
momentum transfer. We see that waves are excited near the 
dip of N 2 (before N 2 rises to a maximum). This is similar 
to the locat ion of wave excitation in real WD models (see 
Section 6.3 1. The dip in iV 2 causes the wave to have a larger 
wavelength in this region, and so it couples to the companion 
star's gravitational potential best in this region of the star. 
Note that although N 2 is smaller near the center of the star, 
no significant wave is excited there since U oc r 2 is negligible. 

Figure [4] shows a plot of F(uj). For this model, F(uj) 
is not a smooth, monotonic function of cj as it is for the 
massive star model studied in Appendix [X] Instead, there 
are many jagged peaks and troughs, causing the value of 



F(uj) to vary by two or three orders of magnitude over very 
small frequency ranges. These features are also present in the 
real WD models, and will be discussed further in Section [7] 
Our numerical results indicate that the peaks of F(oj) can be 
fitted by F(u) oc uj 5 , significantly different from the massive 
star model. 



6.3 Calculation with Realistic White Dwarf 
Model 

We now present our numerical results for tidal excitations in 
realistic WD models. Using the outgoing wave outer bound- 
ary condition, we solved the oscillation equations |5| and |6| 
for the three WD models described in Section [r72| (see Figure 
[T|. Figures [5] and |6] show plots of the outgoing energy flux 
as a function of radius for the model with T e g = 8720K and 
tidal frequencies of uj = 2fl — 10~ 2 and 1.1 x 10~ 2 , in units 
of G = M = R = 1, respectively. The energy flux jumps to 
its final value near the carbon-helium transition zone. Once 
again, the oscillations in energy flux are due to imperfect 
numerical calculation of the dynamical component of the 
wave and do not actually contribute to energy or angular 
momentum transfer. In Figures [5] and [6) we have smoothed 
the value of the energy flux to minimize the amplitude of the 
unphysical oscillations. Note that although Figure [6] corre- 
sponds to a larger tidal frequency, the outgoing energy flux 
is about 100 times less than in Figure[5] Thus, as in our toy 
WD model (see Section 6.2 I, the tidal energy flux is not a 



monotonic function of tidal frequency as it is for early-type 
stars (see Section [A| . 

We have calculated the dimensionless tidal torque 
F(uj) = j z /T [see equation |42| ] as a function of uj for 
the three WD models depicted in Figure [I] The results are 
shown in Figures [7J |8j and [9] In general, F(uj) exhibits 
a strong and complicated dependence on uj, such that a 
small change in uj leads to a very large change in F(uj) (see 
also Figures |5]|6| . This dependence is largely due to "reso- 
nances" between the radial wavelength of the gravity waves 
and the radius of the carbon core, as discussed in Section [7] 
We also find that the local maxima of F(uj) can be approx- 
imately fitted by the scaling F(uj) oc uj , similar to the tow 
WD model discussed in Section [(L2I 



6.4 Relation to Tidal Overlap Integral 

In an attempt to understand the erratic dependence of the 
tidal energy transfer rate E on the tidal frequency uj, here 
we explore the possible relationship between E and the tidal 
overlap integral. The energy transfer rate to the star due to 
tidal interactions can be written as 



E 



-2 Re / d 3 xp d ^ ,t] ■ VU*(r,t). 



dt 



(48) 



With £(r,t) = £(r)e- iwt and U(r,t) 
equations Q-(|2|], we have 



U(r)Y 22 e- iut [see 



• GM'W 22l 
E = 2uj - Im 



d\pi(y)- V(r 2 Y 2 * 2 ; 



(49) 



We decompose the tidal response £(r,i) into the superposi- 
tion of stellar oscillation modes (with each mode labeled by 
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Figure 2. The square of the Brunt Vaisala (thin solid line) and Lamb (dashed line) frequencies (for I = 2), in units of GM/Ft? , as a 
function of normalized radius in a toy WD model. Also plotted is the stellar density (thick solid line). The stellar properties are only 
plotted out to r = 0.9R, the location where the outer boundary condition is imposed in the tidal excitation calculation. 



the index a) 



£(r,t) = X><* (<)€« W, 



(50) 



where the mode eigenfunction £ a is normalized via 
fd 3 x p\£ a \ 2 = 1. Then the mode amplitude a a (t) satisfies 
the equation 



GM'W 22 Q a - iut 
= e 



(51) 



where ui a is the mode frequency, 7^ is the mode (amplitude) 
damping rate, and Q a is the tidal overlap integral with mode 
a: 



Qa = Jd 3 xpC(r)-V(r 2 Y 22 ). 



The steady-state solution of equation (51 1 is 
GM'W 22 Q a 



a a (t) = — 



Thus the tidal energy transfer rate to mode a is 
(GM'W 22 \Q a ^ 1 



E a = 2tv 



J a UJ 



(52) 



(53) 



(54) 



In paper I, we have computed uj a and Q a for adia- 
batic g-modes of several WD models used in this paper. 
The eigenfunctions of these modes satisfy the "reflective" 
boundary condition (i.e., the Lagrangian pressue perturba- 
tion AP vanishes) at the WD surface. Our result showed 



that although the mode frequency uj a decreases as the ra- 
dial mode number n increases (for a given 1 = 2), the overlap 
integral \ Q a \ is a non-monotonic function of n (or cj a ) due to 
various features (associated with carbon-helium and helium- 
hydrogen transitions) in the TV 2 profile of the WD models. 
On the other hand, our calculation of the tidal response 
£(r, t) presented in this paper adopts the radiative outer 
boundary condition; this implies significant wave damping 
at the outer layer of the star. Because of the difference in 
the outer boundary conditions, the mode frequency ui a (as 
computed using the AP = boundary condition) does not 
have special significance. Nevertheless, we may expect that 
when uj — ui a , the tidal energy transfer is dominated by a 
single mode (a) and E is correlated to |Q a | 2 . 

In Figures [T| [8] and [9] we show \Q a \ as a function of 
ui a for a number of low-order g-modes. It is clear that the 
peaks and troughs of F(oj) calculated with an outgoing wave 
outer boundary condition are associated with the peaks and 
troughs in the value of |Q a | 2 . Thus, the peaks in the value of 
F(u>) are not due to resonances with g-modes, but approxi- 
mately correspond to the tidal frequencies near the "intrin- 
sic frequencies" of the g-modes with large values of |Q a | 2 . 
Note this correspondence between |Q Q | 2 and the local peaks 
of F(uu) is not precise (as they are calculated using differ- 
ent boundary conditions), as is clear from the T e g = 2720K 
model (Figure j9|). Another way to understand the erratic 
dependence of F(u) on ui lies in the quasi-resonance cavity 
of the carbon core of the WD (see Section . 




Figure 3. Dynamical tide in a toy WD model (based on the model depicted in Figure [2] driven by a companion of mass M' = M, 
with the t ida l frequency uj = 2.3 X 10 — 2 . Top: The energy flux E = QJ Z (dark solid line) as a function of radius, with J z calculated from 
equation |39| . All values arc plotted in units of G = M = R = 1. Bottom: The real part of £^ yn (dark solid line) and imaginary part of 
£ dyn (dashed line) as a function of stellar radius. The value of N 2 has been plotted (light solid green line) in both panels. In this model, 
the energy flux rises to its final value near the dip in N 2 , showing that the wave is excited at this location. 



6.5 Justification of the Outer Boundary Condition 

Our calculations in this paper adopt the outgoing wave 
boundary condition near the stellar surface. This implicitly 
assumes that gravity waves are absorbed in the outer layer of 
the WD due to nonlinear effects and/or radiative damping. 
To analyze the validity of this assumption, we plot the mag- 
nitude of the displacement, |£ dyn |, as a function of radius in 
Figure |I0| We have shown the results for tidal frequencies 
of uj = 2fl = 0.028 and 0.0053 (corresponding to frequen- 
cies near the peaks in F(ui) shown in Figure [7| for our WD 
model with T e g = 8720K. We have also plotted the local 
radial wavelength k~ x because we expect nonlinear wave 
breaking to occur when |£ dyn | > k~ x . 

It is evident from Figure [10] that at relatively high tidal 
frequencies, the gravity waves become nonlinear in the outer 
layer of the star, justifying our outgoing wave boundary con- 
dition. In some cases, the waves formally reach nonlinear am- 
plitudes {k~ x \^ y ^\ > 1) in the helium-hydrogen transition 
region (demarcated by the dip in k^ 1 at r ~ 0.935). This 
implies that waves may be partially reflected at the helium- 
hydrogen transition region, although nonlinear damping 



may also occur before the waves make it to the outermost 
layers of the WD. The lower frequency gravity waves do not 
formally reach nonlinear amplitudes in the region depicted 
in Figure [10] However, when extending to the lower-density 
region near the stellar surface, the wave amplitudes will in- 
crease further and nonlinearity will set in, although partial 
reflection may occur due to the shallow convection zone very 
near the stellar surface. Also note that our calculations are 
for Q„ = 0. If the WD has a non-negligible spin (fi s ), a given 
tidal frequency uj — 2(fi — Q s ) would correspond to a higher 
orbital frequency Q, further increasing the wave amplitudes 
compared to those shown in Figure |l0| Furthermore, lower 
frequency waves may damp efficiently via radiative diffusion 
near the stellar surface. We therefore expect our outgoing 
wave outer boundary condition to be a good approximation 
for the frequencies considered in this paper for our warmest 
WD model. 

Our cooler WD models with T cH = 5080K and T cff = 
2720K do not formally reach the same nonlinear ampli- 
tudes as our warmest model. The cooler models have smaller 
Brunt- Vaisala frequencies, particularly in their outer lay- 
ers, as can be seen in Figure [I] Consequently, the gravity 




0.01 0.10 

Figure 4. The dimensionless tidal torque F(ui) = J z /T [see equation fl42J|] carried by the outgoing gravity wave as a function of the 
tidal frequency u) (solid line), for the toy WD mod el d epicted in Figure [2] The frequency is in units of G = M = R = 1. The straight 
light solid (green) line is calculated from equation j68| l and is roughly proportional to ui 5 . The dashed (red) line is our semi-analytical 
approximation, with a = 1/5, /3 = 1/5, and 5 = (see Section[7]l. 



waves have smaller displacements (recall the WKB scaling 
£^ yn oc iV 1,/2 for a constant j z ) and larger wavelengths (re- 
call k r oc iV). Therefore, gravity waves are less likely to 
damp due to nonlinear effects in our cooler models, and our 
outgoing wave outer boundary condition may not be justi- 
fied at all frequencies considered. More detailed analyses of 
the nonlinear effects in dynamical tides are necessary (e.g., 
Barker & Ogilvie 2010, Weinberg et al. 2011). 



with 



k 2 (r) 



Z" + k 2 (r)Z = V(r), 
l(l + l)N 2 



l(l + l)N 2 



+ Ar(r) 



1 + 



H 2 l(l + l)N 2 



(55) 



(56) 



7 SIMPLE MODEL FOR GRAVITY WAVE 
EXCITATION: ANALYTICAL ESTIMATE 

To understand our num erica l result for the tidal energy 
transfer rate E (Section 6.3 1, particularly its dependence 



on the tidal frequency co, here we consider a simple stellar 
model that, we believe, captures the essential physics of tidal 
excitation of gravity waves in binary WDs. In this model, 



the star consists of two regions (see Figure 111: the outer re- 
gion with r > r a (region a) and the inner region with r < r;, 
(region b). In each region, the stellar profiles are smooth, 
but TV 2 jumps from N 2 at r = r b to N 2 (with TV 2 > N 2 ) 
at r — r a . The tidal frequency u> satisfies u> 2 <C iV 2 . As we 
will see, although waves can propagate in both regions, the 
sharp jump in iV 2 makes the inner region behave like a res- 
onance cavity-this is ultimately responsible for the erratic 
dependence of F(co) on the tidal frequency ui. 



We start from the wave equation ( 12 1 for Z(r) 



V{r) 



_ - 1/2 l{l+l)N 2 U 
to 2 g 



2r 



H l{l + l)N 2 



(57) 



where H = a 2 /g(& r) is the pressure scale height. The above 
expressions are valid in both regions of the star, and we have 
assumed to 2 <g L 2 and uj 2 <^ N 2 [more general expressions 



are given by equations (13 1 and ( 14 1] . Since the stellar pro- 



files are smooth in each of the two regions, the non-wave 
("equilibrium") solution is given by 



where 



<7 -1/2 



r 2 U 

1 

9 



AZ = Zo 



k 2 H 2 



(58) 

(59) 
(60) 
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Figure 5. Dynamical tide in a realistic WD model (with M = 0.6M Q , R = 8.97 X 10 3 km, and T off = 8720K) driven by a companion 
of mass M' = M, with the tidal frequency u = 2Q = 10~ 2 . Top: The energy flux E = QJ Z (thick solid line) as a function of radius, 
calculated from equation |39| , All values are plotted in units ofG = M = iJ=l. Middle: The real part of £^ yn (solid line) and imaginary 
part of £^ yn (dashed line) as a function of stellar radius. Bottom: The same as the middle panel, but zoomed in on the outer layer of 
the WD. The value of N 2 has been plotted as dashed (green) lines in each panel. The energy flux rises to near its final value around the 
carbon-helium transition region, showing that the wave is excited at this location. 



and (3 is a constant (with /3 ~ 1)|^] Note that the above 
solution for Z eq breaks down around r = r- m (where lu 2 — 



3 From equation ( |13| 
byAk 2 = (H~ 1 Y/2- 
where H p = H and 1 
H, and we have Ak 2 



p 5 / 3 , we have H p = (5/3)H, and Ak 2 
parameter /3 in equation i59\ ranges from 



, we find that Ak 2 in equation (1561 is given 
■ (2H P )~ 2 + H^l-ilnHp)' + HJ* ~ H' 1 }, 
—p'/p. In the isothermal region, H p = 
-(2H P )~ 2 . In the region satisfying P oc 



l/(12_ff2). Thus the 
31 < 0.1 to 0.3. 



A^ 2 ). At distances sufficiently far away from ri n , we have 
fc>l/#. 

The general solution to equation \12\ consists of the 
non-wave part Z eq and the wave part Z dyn . In region b there 
exist both ingoing and outgoing waves. Thus 



Z(r) =Z ccl (r) + A + exp^ij kdr 
+ A- exp I —i I kdr 



(61) 
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Figure 6. Same as FigureJE] but for the tidal frequency w = 1.1 X 10 



for ri n < r < rb, where A+ and A- are slow- varying func- 
tions of r. In region a, we require there be no ingoing wave. 
Thus for r > r a , 



Z(r) = Z cq (r) + Aexp 



il kdr 



(62) 



where ro > r a is a constant, and A varies slowly with r. 
Note that Z eq is discontinuous between the two regions. 

At the inner boundary r = r- m , gravity waves are per- 
fectly reflected. Thus we have A- = —e lS A+, where S is 
a constant phase that depends on the details of the dis- 
turbance around and inside n n . To determine A and A+ 
we must match the solutions in the two regions. Although 
in reality r a is somewhat larger than rt, we shall make 
the approximation r a — r~b, and label the physical quan- 



tities on each side with the subscript "a" or "b". Note 
that - Z cq ~ -AZ b since k 2 a > fcg, and (dZ ccl /dr) a - 
(dZ a<1 /dr) b ~ —(a/H)AZ b where a is a constant (|a| ~ 1). 
Matching Z and dZ/dr across r — r b — r a , we obtain the 
expression for the wave amplitude at r — r a : 



A exp 
which entails 
\A 

where 



z I kdr 

ra 



AZ, 



(a/k b H) tan <p 



\AZ b 



_ 1 + i{k a /k b )\axitp 
(a/kbH) tan ip\ 



[l + (fc a /fc b ) 2 tan 2 v9] 1/2: 



<P = 



kdr — - . 
2 



(63) 



(64) 



(65) 
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Figure 7. The dimensionless tidal torque F(u) = j z /T [see equation (42| ] carried by outgoing gravity waves as a function of the 
tidal frequency w for our WD model with T e g = 8720K. The two boxed points correspond to w = 10 
Figures [5] and [6] The dashed (red) line is our semi-analytical approximation [see equation 167} ], with a 
smooth solid line corresponds to the maximum values of F{w) in our semi-analytical equation, and is calculated from equation (69) . The 
dot-dashed (blue) line corresponds to F(ui) = 20(1> 5 (see Section [sj. The diamonds connected by the dotted line are the tidal overlap 
integrals Q a associated with nearby gravity modes, and the n = 4 mode is the highest frequency mode shown. The frequency and Q a 
are plotted in units oiG = M = R= l. 



2 and 1.1 X 10 2 , as depicted in 
1/5, p = 1/5, and 5 = 0. The 



Clearly, \A\ reaches the maximum |AZt| at ip = 0, and \A\ ~ 
\AZb(a/k a H)\ at <p = tt/2. 

The Lagrangian displacement for the outgoing gravity 
wave in region a is given by 

/■dyn ikT j-dyn 

ik X 1/2 



-Aexp 



-I I kdr 

ro 



(66) 



l(l + l)r 

The tidal energy transfer rate E is equal to the energy flux 



carried by the wave. Using equation (41 1, we have 



E = QJ Z = 4Qk a \A\ 2 , (67) 

where k a = \fl{l + l)N a / {r a uj) and \A\ are evaluated at 
r — r a . Using |A| max = |A.Z&|, we obtain the maximum tidal 
energy transfer rate as a function of the tidal frequency lj 
and the orbital frequency fl: 



E n 



37T/3 2 



ParlNa 



N£ll(l + l)] 5/2 g%\H a 



\ f M' \ o5 .5 



(68) 



The corresponding dimensionless tidal torque [see equation 
(42)] is 

Gp a r 7 a N a 



F n - 



6tt/3 2 



5 Nt[l(l + l)f /2 gl& \H a 



(69) 



This scaling [F(u) oc ui 5 ] agrees with our numerical results 



for the toy WD models (Section 6.2 1 and realistic WD mod- 
els (Section 



■•3i 



Realistic WD models are obviously more complicated 
than the analytical model considered in this section (see 
Figure 11 1 . To evaluate the tidal energy trans fer rate E using 



equation (67 1 [with \A\ given by equation ( 64 1] and E-, 



using equation (68 1 for our WD models, we choose n, at the 
location where d\nN 2 /dr is largest in the helium-carbon 
transition region. We then set the location of r a to be one 
half of a wavelength above rt, i.e., by finding the location r a 
such that the equation 7t = kdr is satisfied, where k is 



given by equation (271. For the three models considered in 



Section 6, we find that r a thus calculated typically lies near 
the peak in iV 2 associated with the carbon-helium transition 
region. 

In Figures |4j [7] |U andj9l we compare the analytical 
results based on equations (67) and (68) to our numerical 
calculations. We see that the erratic dependence oi F(uj) on 
the tidal frequency lu can be qualitatively reproduced by our 
analytical expression (67), and the maximum F maj( is also 



well approximated by equation ( 69 1. Our analytical estimate 
works best for the WD model with T Q g = 8720K, but it does 
a poor job of approximating the value of F(uj) for the WD 
model with T e ff — 2720K. We attribute this disagreement 
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Figure 8. Same as Figure[7] except for the T e g = 5080K WD model. In this plot, the dot-dashed (blue) line corresponds to F(uj) = 200(i 5 . 
The n = 15 mode is the highest frequency g-mode shown. 



to the lower value of JV in the cool WD model because 
our assumption that iV 2 3> is not satisfied. Instead, we 
find that gravity waves are excited near the spike in N 2 
associated with the helium-hydrogen transition region in the 
cool WD model. 

For each model shown in Figures |7|9| our model also 
breaks down at the highest and lowest frequencies shown. 
These discrepancies are likely related to errors in our nu- 
merical methods. At the highest frequencies shown, the ap- 
proximation k r 3> 1/H begins to break down, causing error 
in our outer boundary condition. At the lowest frequencies 
shown, extremely fine grid resolution is needed to resolve 
the dynamical component of the tidal response, and so slight 
thermodynamic inconsistencies may introduce significant er- 



rors (see Section 6.1l 



8 SPIN-ORBIT EVOLUTION 

The tidally-excited gravity waves and their dissipations 
cause energy and angular momentum transfer from the orbit 
to the star, leading to spin-up of the WD over time. In this 
section, we study the spin-orbit evolution of WD binaries 
under the combined effects of tidal dissipation and gravita- 
tional radiation. In general, the tidal torque on the primary 
star M from the companion M' and the tidal energy transfer 
rate can be written as [see equations (|42[) and (|43[)] 



(70) 



Ttido = T F(oj), iStidc = T QF(w), 



with To = G(M' /a 3 ) 2 R 5 . In previous sections, we have com- 
puted F(uj) for various non-rotating (Q s = 0) WD models 
(and other stellar models). To study the spin-orbit evolu- 
tion, here we assume that for spinning WDs, the function 
F(lj) is the same as in the non-rotating case. This is an ap- 
proximation because a finite Q s can modify gravity waves 
in the star through the Coriolis force (gravity waves become 
the so-called Hough waves) and introduce inertial waves, 
which may play a role in the dynamical tides. In other 
words, the function F generally depends on not only uj but 
also Q s . However, we expect that when the tidal frequency 
uj — 2(Q — Q s ) is larger than Q s , i.e., when Q <; 3Q 3 /2, the 
effect of rotation on the gravity waves is small or modest. 
Also, we assume that the WD exhibits solid-body rotation, 
which would occur if different layers of the WD are strongly 
coupled (e.g., due to viscous or magnetic stresses)]^] 

Before proceeding, we note that in the weak friction 
theory of equilibrium tides (e.g., Darwin 1879; Goldreich & 
Soter 1966; Alexander 1973; Hut 1981), the tidal torque is 
related to the tidal lag angle <5i ag or the tidal lag time A£i ag 



4 In a medium containing a magnetic field, we expect differen- 
tial rotation to be smoothed out by magnetic stresses on time 
scales comparable to the Alfven wave crossing time. The Alfven 
wave crossing time is £a = R\J4>np/B ~ lyr for a magnetic field 
strength of B = 10 5 gauss and a density of p = 10 6 g/cm 3 . Since 
the Alfven wave crossing time is always much smaller than the 
inspiral time for WDs, we expect solid body rotation to be a good 
approximation. 
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Figure 9. Same as Figure [7] except for the T e g = 2720K WD model. In this plot, the dot-dashed (blue) line corresponds to F(ui) 
4 X 10 3 (2> 5 . The n = 9 mode is the highest frequency g-mode shown. 



by 



Ttide = 3feT (5i a g, with 5i ag = (fi — fi s )Ai- 



lag! 



(71) 



where fc 2 is the Love number. Often, a dimensionless tidal 
quality factor Qtidc is introduced (e.g. Goldreich & Soter 
1966) such that Aii ag = l/(|w|<3tide) (valid only for to 0). 
Thus, if we use the weak-friction theory to parametrize our 
dynamical tide, F(lj) would correspond to 

F(uj) = 3fc 2 (5lag 

= 3fc 2 (^-« 3 )Ati ag 
3k 2 



sgn(fi — il s 



(72) 



Obviously, the effective Qtido would depend strongly on ui as 
opposed to being a constant (assuming constant lag angle) 
or being proportional to l/|w| (assuming constant lag time, 
appropriate for a viscous fluid). 

With equation (701 and the assumption in F(uj), the 



WD spin evolves according to the equation 



(73) 



where I is the moment of inertia of the WD (7 ~ 0.169M R 2 
for our M = 0.6M Q WD models). The orbital energy E orb = 
—GMM'/(2a) satisfies the equation 



-Etide — Eg\ 



(74) 



where Egw(> 0) is the energy loss rate due to gravitational 
radiation. The evolution equation for the orbital angular 



frequency Q = {GM t /a 3 ) 1/2 is then 



3T F(u) 30. 



2tGW 



(75) 



where \i — MM' /Mt is the reduced mass of the binary, and 
taw is the orbital decay time scale (|a/d|) due to gravita- 
tional radiation: 



5c : ' 



64G 3 MM'M t 

3.2 x lO 1 ^ M ^ ( M < V 



\MM'J\2MqJ \0.1 



-8/3 



s, (76) 



8.1 Synchronization 

Using our results for the function F(lj) obtained in previous 



sections, we integrate equations ( 73 1 and ( 75 1 numerically to 



obtain the evolution of the WD spin. Since at large a (small 
57) the orbital decay time ~ taw oc f2~ 8//3 is is much shorter 
than the time scale for spin evolution, t sp in = O s /Q s oc 
1/(Q 4 F), we start our integration with fl s < O at a small 
orbital frequency (an orbital period of several hours). 

The results for our three WD models are shown in Fig- 
ures [12] and [13] Note we only include the effects of tides 
in the primary star (M), and treat the companion (M') as 
a point mass. All three models have the same WD masses 
(M = M' = O.6M0), but different temperatures. Also note 
that the minimum binary separation (before mass transfer 
or tidal disruption occurs) is a m i n ~ 2.5(Mt/M) 1 ^ 3 R, corre- 
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Figure 10. The magnitude of the gravity wave displacement vector |£ d y n | (solid line) as a function of radius for a tidal frequency of 
u = 0.028 (top) and w = 0.0053 (bottom). Also plotted is the local radial wavelength (red dashed line). The wave displacement, 
wavelength, and frequency are in units where G = M = R = 1. 



sponding to the minimum orbital period 



P min ~ (1.1mm) M~ 1/2 R 3 4 /2 , 



(77) 

where Mi = M/(1M ) and R 4 = 7?/(10 4 km). We see 
that for all models, appreciable spin-orbit synchronization is 
achieved before mass transfer or tidal disruption. However, 
depending in the WD temperature, the rates of spin-orbit 
synchronization are different. 

The basic feature of the synchronization process can be 
obtained using an approximate expression for the dimension- 
less function F(u>). We fit the local maxima of our numerical 
results depicted in Figures [7|9| by the function 

(78) 



where u = ui/(GM/R 3 ) 1/2 , and / ~ 20, 200, 4 x 10 3 for 
the T cS =8720K, 5080K, and 2720K models, respectively. 
Suppose Q. s <g; Q. at large orbital separation. We can de- 
fine the critical orbital frequency, £l c , at which spinup or 
synchronization becomes efficient, by equating Cl and Cl s 
(with f2 s <C £2). Note that since the orbital decay rate due 



to tidal energy transfer [the first term in equation ( 75 1] 
is much smaller than the spinup rate Cl s , the orbital de- 
cay is always dominated by the gravitational radiation, i.e., 
f2 ~ 3fi/(2*Gw). With T = f Q, A and t GW = t GW D.- 8/3 , 



F(w) = fuj 5 = fuj 5 , 
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Figure 11. A diagram showing a simplified model of a white dwarf used in our analytical estimate for gravity wave excitation. The arrows 
indicate that region b contains both an inward and outward propagating wave, while region a contains only an outward propagating 
wave. 



we find 



3/ 



64/TUgw 

r 5/3 



3/16 
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5/2' 



3/16 



(GM\ 1/2 
\B? ) 



(3.8 x 10~ 3 s _1 ) 
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wo.i 7 M t Y a M 1 a 
fM[Rf /2 



3/16 



(79) 



where n = 0.17« .i7 = I /{MR 2 ), M{ = M'/(1M Q ), and 
M t i = Mt/(1 Mq). For Q ^ £} c , tidal synchronization is in- 
efficient. For f2 <; fi c , the spinup rate Cl 3 becomes larger than 
Q and the system will become increasingly synchronized. In 
fact, when Q <; £l c , an approximate analytic expression for 
the spin evolution can be obtained by assuming a posteriori 
(fi a - fi) < f2. With f2 ~ 3fi/(2t G w) ^ we find 



16/15^-1/15 



(for n ;> n c ). 



(80) 



This expression provides an accurate representation of the 
numerical solutions. 



Note that we can derive a similar equation as ( 80 1 for 
more general tidal torques. For example, assume 

Qs = AO." 1 ^ -Q s ) n , (81) 

where n and A are constants. With £l = BQ 11 ^ 3 (where B 



is a constant) and assuming Q s ~ £1, we find 



l/(3n) 



for f2 <; f2 c , where 



a. 



(3) 



3/(3n+l) 



(82) 



(83) 



Note that our equation ( |80[ ) corresponds to n = 5, which 
implies £1 — f2 a ~ f2 c for fi <; fi c . By contrast, in the equilib- 
rium tide model (with constant lag time), n = 1, so (fi — Q s ) 
changes moderately as the orbit decays. 



8.2 Tidal Effect on the Orbital Decay Rate and 
Phase of Gravitational Waves 



Figure 14 shows the tidal energy transfer rate (from the 
orbit to the WD) £Wc = T QF(u). For Q £ fi c , cj ~ 2Q, 
(assuming Q. s <C Q), we see that -Etidc depends on f! in a 
rather erratic manner. However, when fi <; f2 c , efficient tidal 
synchronization ensures f2 ~ fi s , or 3f2/2iGW — ToF(u)/I, 
and thus Stide simplifies to 



-Etid 



Since E tide /E c 



2tow 
3//(/V) 



(forfi£fi c ). 



(84) 



<C 1, the orbital decay is dom- 



inated by gravitational radiation. Nevertheless, the orbital 
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Figure 12. Evolution of the spin frequency Q s in units of the orbital frequency Q as a function of the orbital period. The solid black, 
red, and blue lines correspond to our WD models with T e g = 8720K, T = 5080K, and T = 2720K, respectively. The black, red, and 
blue dashed lines correspond to evolutions using F = 20a; 5 , F = 200lD 5 , F = 4x 10 3 aj 5 , respectively (these functions F(oj) approximate 
the like-colored WD models, see Figures |7)9| l. The vertical dotted line denotes the critical orbital period, 2tt/Q c [see equation | |79| l], 
corresponding to the black dashed line. In these evolutions, M' = M and the WDs initially have Q s = 0. 



phase evolution is affected by the tidal energy transfer, and 
such a phase shift can be measurable for short period bina- 
ries such as the recently discovered 12 minute system SDSS 
J0651 (Brown et al. 2011; see Section 9). Also, low-frequency 
(10 -4 — 10 _1 Hz) gravitational waveforms emitted by the bi- 
nary, detectable by LISA, will deviate significantly from the 
point-mass binary prediction. This is in contrast to the case 
of neutron star binaries (NS/NS or NS/BH) studied pre- 
viously (Reisenegger & Goldreich 1994; Lai 1994; Shibata 
1994; Ho & Lai 1999; Lai & Wu 2006; Flanagan & Racine 
2007), where the resonant mode amplitude is normally too 
small to affect the gravitational waveforms to be detected 
by ground-based gravitational wave detectors such as LIGO 
and VIRGO, tidal effects only become important near the 
NS binary merger (e.g., Lai et al. 1994; Hinderer et al. 2010). 
The orbital cycle of a WD binary evolves according to 

dA r orb = _p_^. (85 ) 
2^ E OTh 



Including tidal effects in i5 or b, we find 



diVorb / diVorb 



dlnfi \dlnf2/ Q \ Eg 



1 + 



Etid 



(86) 



where 

( dN OTb 



r>c" 



\dlnftj 3tt 192TrG 5 / 3 nM? /3 (nf G w) 5/3 



2.3x10 



/ A/f2 \ / nr \V 3 / r \-S/ 3 
9 / M Q \ I M t \ ( JGW N 



V MM' J V 2Mq J \0-01Hz 



(87) 



is the usual result when the tidal effect is neglected (/gw = 
Q/n is the gravitational wave frequency). Thus, even though 
Etide/EGw <§C 1, the number of "missing cycles" due to the 
tidal effect, 



dAN olb 
dlnft 



dN OT b \ -Etidc 



d In Q J _E GW 



(88) 



can be significant. Since -E t id G oc I, proper modelling and 
detection of the missing cycles would provide a measurement 
of the moment of inertia of the WD. 



8.3 Tidal Heating 

The tidal energy transfer -Etide does not correspond to the 
energy dissipated as heat in the WD, because some of the 
energy must be used to spin up the WD. Assuming rigid- 
body rotation, the tidal heating rate is 



-Eheat — -Etid 



1 - 



(89) 



Figure 



14 



shows E?heat for our three binary WD models. At 
large binary seperations (Q <i Q, c ) when £l s <C fi, virtu- 
ally all of the tidal energy transfer to the WD is dissipated 
as heat. At smaller serparations, we have shown that the 
WD will retain a small degree of asynchronization. Insert- 
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Figure 13. The spin frequency Q s /(2tt) in units of Hz as a function of orbital period. The solid black, red, and blue lines correspond 
to our WD models with T cff = 8720K, T = 5080K, and T = 2720K, respectively. The dashed line shows the orbital frequency, U/(2tt). 
The dotted vertical black, red, and blue lines are the values of 27r/f2 c for F = 2(K> 5 , F = 20(K> 5 , F = 4x 10 3 £ 5 , respectively. In these 
evolutions, M' = M and the WDs initially have U s = J7/4. 



ing equation (1801) into equation i89\, we find 



9 DISCUSSION 



E h 



-Etid 



3m 2 fo. 



2^gw V ^ 



16/15 



16/15 



(for n z n c 



(90) 



Thus, as the orbital frequency increases, a smaller fraction 
of the tidal energy is dissipated as heat. Using equation ( 79 1 
for fi c , we have 



T? fc i , io36 -1\ 6/5 ?-l/5 , j-29/10 

Sheat ~(6.1 x 10 erg s )k . 17 / ' M x 
/ o \ 18 / 5 

x (Af0 4/5 ^ 4 " 1/10 



0.1 s" 



(91) 



Note that Eheat is relatively insensitive to /, so its pre- 
cise value is not important. Thus, tidal heating of the WD 
can become significant well before merger. For example, for 
our Tetf = 8720K WD model (with M = M' = 0.6M o , 
R = 8970 km and / ~ 20), we find E he3lt ~ 1.2 x 10 32 erg/s at 
the orbital period P = 10 min, much larger than the "intrin- 
sic" luminosity of the WD, 47ri? 2 CT S BT , c 4 ff = 3.3 x 10 30 erg/s. 
Note that E^cat is mainly deposited in the WD envelope, so 
an appreciable fraction of £"tide may be radiated, and the 
WD can become very bright prior to merger. The 12 minute 
binary SDSS J0651 (Brown et al. 2011) may be an example 
of such tidally heated WDs (see Section 9). 



We have studied the tidal excitation of gravity waves in bi- 
nary white dwarfs (WDs) and computed the energy and an- 
gular momentum transfer rates as a function of the orbital 
frequency for several WD models. Such dynamical tides play 
the dominant role in spinning up the WD as the binary 
decays due to gravitational radiation. Our calculations are 
based on the outgoing wave boundary condition, which im- 
plicitly assumes that the tidally excited gravity waves are 
damped by nonlinear effects or radiative diffusion as they 
propagate towards the WD surface. Unlike dynamical tides 
in early-type main-sequence stars, where gravity waves are 
excited at the boundary between the convective core and 
radiative envelope, the excitation of gravity waves in WDs 
is more complicated due to the various sharp features as- 
sociated with composition changes in the WD model. We 
find that the tidal energy transfer rate (from the orbit to 
the WD) .Etidc is a complex function of the tidal frequency 
uj — 2(0 — Q s ) (where fi and Q s are the orbital frequency 
and spin frequency, respectively; see Figures 7]9 1, and the lo- 



cal maxima of E t id B scale approximately as fl or . For most 
tidal frequencies considered, the gravity waves are excited 
near the boundary between the carbon-oxygen core and 
the helium layer (with the associated dip and sharp rise in 
the Brunt- Vaisala frequency profile). We have constructed a 
semi-analytic model that captures the basic physics of grav- 
ity wave excitation and reveals that the complex behavior of 
E t idc £ls a function of the tidal frequency arises from the par- 
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Figure 14. The tidal energy dissipation rate -Etide (solid lines) and the tidal heating rate -Ehcat (dashed lines) as a function of orbital 
period. The black, red, and blue lines correspond to our WD models with T e g = 8720K, T = 5080K, and T = 2720K, respectively. Note 
that at small orbital periods, the i?tidc curves overlap for different WD models. The dotted line is the energy dissipation rate due to 
gravitational waves, Egw- In these evolutions, M' = M and the WDs initially have Q s = 0. 



tial trapping of gravity waves in the quasi-resonance cavity 
provided by the carbon-oxygen core. 

We have also calculated the spin and orbital evolution 
of the WD binary system including the effects of both gravi- 
tational radiation and tidal dissipation. We find that above a 
critical orbital frequency Q c [see equation ( 79 1], correspond- 



ing to an orbital period of about an hour for our WD mod- 
els, the dynamical tide BEGINS to drive the WD spin Q s 
towards synchronous rotation, although a small degree of 
asynchronization is maintained even a t sm all orbital periods: 
Q — O s ~ fi c (fi c /f2) 1,/15 [see equation ( 80 1] . Thus, numerical 



simulations of WD binary mergers should use synchronized 
configurations as their initial condition - these may affect 
the property of the merger product and possible supernova 
signatures. 

We also show that, although gravitational radiation al- 
ways dominates over tides in the decay of the binary orbit, 
tidal effects can nevertheless affect the orbital decay and 
introduce significant phase error to the low-frequency gravi- 
tational waveforms. Future detection of gravitational waves 
from WD binaries by LISA may need to take these tidal 
effects into account and may lead to measurements of the 
WDs' moments of inertia. Finally, we have calculated the 
tidal heating rate of the WD as a function of the orbital pe- 
riod. For Q <; Sl c , since the tidal dissipation rate is largely 
controlled by the orbital decay rate due to gravitational ra- 
diation, it is a smooth function of orbital period [see equa- 



tion (90 1], We show that well before mass transfer or binary 
merger occurs, tidal dissipation in the WD envelope can be 
much larger than the intrinsic luminosity of the star. Thus, 



the WD envelope may be heated up significantly, leading to 
brightening of the WD binary well before merger. We plan 
to study this issue in detail in a future paper. 

The recently discovered 12 minute WD binary SDSS 
J0651 (Brown et al. 2011) can provide useful constraints for 
our theory. Applying equation ( 79 1 to this system, we find 



that the orbital period (12.75 minutes) is sufficiently short 
that both WDs are nearly (but not completely) synchronized 
with the orbit. Because of the orbital decay P, the eclipse 
timing changes according to the relation 



At = Pt 2 /{2P), 



(92) 



where t is the observing time. Gravitational r adia tion gives 
rise to /S±gw = 5.6s (i/lyr) 2 . Using equation (841 to evalu- 
ate the orbital decay rate Ptide due to tidal energy transfer, 
we find A£ t id c — 0.28s (i/lyr) 2 (see also Benacquista 2011). 
Thus, the orbital decay due to tidal effects should be mea- 
sured in the near future. Also, our calculated heating rate, 
equation | |91| |, indicates the SDSS J0651 WDs have suffered 
significant tidal heating, although to predict the luminos- 
ity change due to tidal heating requires careful study of the 
thermal structure of the WDs and knowledge of the location 
of tidal heating. We note that Piro (2011) also considered 
some aspects of tidal effects in SDSS J0651, but his results 
were based on parameterized equilibrium tide theory. 

This paper, together with paper I, represents only the 
first study of the physics of dynamical tides in compact WD 
binaries, and more works are needed. We have adopted sev- 
eral approximations that may limit the applicability of our 
results. First, we have not included the effects of rotation 
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(e.g., the Coriolis force) in our wave equations. In addition 
to modifying the properties of gravity waves (they become 
generalized Hough waves), rotation also introduces inertial 
waves that can be excited once the WD spin frequency be- 
comes comparable to the tidal frequency - this may lead to 
more efficient tidal energy transfer and synchronization. For 
example, if we parameterize the spinup rate due to various 
mechanisms (including inertial waves) by equation (87), the 
critical orbital frequency for the onset of synchronization 
(f2 c ) is given by equation (831. For a stronger tidal torque 



(larger A), f2 c is smaller. However, the tidal heating rate at 
Q Jj fl c becomes [cf. equation ( 90 1] 



tGW 



(3n+l)/(3n) 



(93) 



Thus, for stronger tidal torques, at a given orbital frequency 
(Q <; fl c ), the tidal heating rate is reduced because the WD 
is closer to synchronization. 

Second, we have assumed that the WD rotates as a rigid 
body. As the tidally-excited gravity waves deposit angular 
momentum in the outer layer of the WD, differential rota- 
tion will develop if the different regions of the WD are not 
well coupled. Thus it may be that the outer layer becomes 
synchronized with the companion while the core rotates at 
a sub-synchronous rate, analogous to tidal synchronization 
in early-type main-sequence stars (Goldreich & Nicholson 
1989). Third, we have implicitly assumed that the outgoing 
gravity waves are efficiently damped near the WD surface. 
This may not apply for all WD models or all orbital frequen- 
cies. If partial wave reflection occurs, tidal dissipation will 
be reduced compared to the results presented in this paper 
except when the tidal frequency matches the intrinsic fre- 
quency of a g-mode (cf. Paper I). More detailed studies on 
nonlinear wave damping (e.g., Barker & Ogilvie 2010; Wein- 
berg et al. 2011) and radiative damping would be desirable. 

Finally, we have only studied carbon-oxygen WDs in 
this paper. Our calculations have shown that the strength 
of dynamical tides depends sensitively on the detailed in- 
ternal structure of the WD. Recent observations (see ref- 
erences in Section 1) have revealed many compact WD bi- 
naries that contain at least one low-mass helium-core WD. 
The temperatures of these helium-core WDs tend to be high 
(T e ff <; 10 4 K). These observations warrant investigation of 
tidal effects in hot, helium-core WDs, which have signifi- 
cantly different internal structures from the cool, carbon- 
oxygen WDs considered in this paper. 
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APPENDIX A: CALCULATION WITH 
MASSIVE STAR MODEL 

To test the accuracy of our numerical calculations and espe- 
cially the importance of self-consistency in real stellar mod- 
els, we compute the tidal response of several toy models. The 
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Figure Al. The square of the Brunt- Vaisala (thin solid line) and Lamb (dashed line) frequencies (for I = 2), in units of GM/R 3 , as a 
function of the normalized radius in a simple massive star model. Also plotted is the stellar density profile (thick solid line) in units of 
M/R 3 . The model has an inner convection zone extending to r = 0.25R. The stellar properties are only plotted out to r = 0.6R, where 
an outgoing wave boundary condition is adopted in our calculation of the tidal excitation. 



first toy model we employ is shown in Figure |AT] and is meant 
to mimic a massive early-type star. The model contains an 
inner convection zone surrounded by a thick radiative enve- 
lope. The convection zone extends to r = 0.25R, beyond 
which the value of N rises linearly to N 2 « 8GM/R 3 . 
Dynamical tides in such massive stars have been studied 
by Zahn (1975, 1977) and Goldreich & Nicholson (1989), 
who showed that the dominant effect arises from the grav- 
ity waves launched at the core-envelope boundary, which 
then propagate outwards and eventually dissipate near the 
stellar surface. Zahn (1975) derived an analytic solution for 
the wave amplitude and the corresponding tidal torque. Al- 
though our model does not contain some of the details ex- 
hibited by realistic massive star models, it does capture the 
most important features. We can compare our result with 
Zahn's to calibrate our numerical method and to assess the 
degree of self-consistency required to produce reliable results 
for the tidal torque. 

Figure |A2| shows an example of our numerical results 
for the dynamical tides generated in a massive star by a 
companion, for a given tidal frequency to — 2Q = 2.3 x 1CP 2 
(in units where G = M — R). We see that gravity waves 
are excited at the base of the radiative zone where N 2 be- 
gins to rise above zero. A net energy flux E = Q.J Z — 
Q(GM' 2 R 5 /a 6 )F(ui) flows outwards toward the stellar sur- 
face. Figure [A3] shows our numerical result of the dimension- 
less function F(lj) = J z /T [see equation (42 1], evaluated at 



,8/3 



in agreement with 



The result can be fitted by F(ui) oc uj~ 
the scaling found by Zahn (1975). 

The power-law scaling of the energy flux can be derived 
using the method of Goldreich & Nicholson (1989). Assume 
|£j_ yn | « at r = r c +, which is located one wavelength 
above the convective boundary (r = r c ). From the dispersion 
relation (271, we find that the Brunt- Vaisala frequency at 



c + is given by (for I = 2) 



N(r c 



dN 2 
dr 



1/3 



M 3 



(Al) 



Using ^ q 
equation (41| to find 



•[l/(6r)] {Ur 2 /g)' ~ -U/(2g), we evaluate 



Ei 



3ttV6 (M' 
10 \Mt 



8/3 



1 pr 7 n 5 oj> 
g 2 (dN 2 /dlnry/ 3 ' 



(A2) 



the outer boundary, as a function of the tidal frequency u>. 



where Mt = M + M' , and all the quantities (p, r, g, and 
dN 2 /dr are evaluated at r = r c + ~ r c ). The scaling of this 
estimate nearly agrees Goldreich & Nicholson (1989), who 
obtained E r oc Sl 4 ^; 11 ^ 3 , where E r is the energy flux carried 
by outgoing gravity waves in the rotating frame of the star, 
not the total energy transfer rate from the orbit. These two 
energy transfer rates are related by E — QJ Z — 2QE r /uj. 
Goldreich & Nicholson (1989) estimates dN 2 /dr w g/H w 
g/r, and with g ~ 4nGpr/3 (p is the mean density interior 
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Figure A2. Dynamical tide in a massive star (based on the toy model depicted in Figure [Al} driven by a companion of mass M' = M, 
with the t ida l frequency u) = 2.3 X 10 — 2 . Top: The energy flux (dark solid line) E = QJ Z as a function of radius, with J z calculated from 
equation |39| . All values are plotted in units of G = M = R = 1. Bottom: The real part of £^ yn (dark solid line) and imaginary part 
of S,^ 11 (dashed line) as a function of stellar radius. The value of N has been plotted in green (light solid line) in both panels. In this 
model, the energy flux rises to its final value just outside of the convective zone, showing that the wave is excited at this location. 



to r c ), equation (IA2I) becomes 



The value of F(ui) based on equation ( |A2[ ) is plotted in Fig- 
ure |A3| Compared to our numerical results, we see that 



equation (A2| overestimates F(io) by an order of magni- 
tude [by contrast, equation ( |A3[ ) would overestimate F(uj) 
by more, since for our toy stellar model dN 2 /dlnr 2> g/r]. 
From our numerical results, we find that the dynamical part 
of the tide only reaches an amplitude of £^ yn w £l q /4- If we 



had used this wave amplitude in our estimate, equation ( A2 1 
would be a factor of 16 smaller and would provide an accu- 
rate approximation to F(oj) at all frequencies considered. 
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Figure A3. The dimensionless tidal torque F(uj) = J z /T [see equation ([42k] carried by the outgoing gravity wave as a function of 
the tidal frequency oj (solid line). The analytical estimate from equation ( |A2| is also plotted (dashed line). The frequency is in units 
ofG = M = i?=l. The small wiggles at high frequencies are likely due to the slight inaccuracy of our implementation of the outer 
boundary condition due to the neglected terms which become non- negligible at higher tidal frequencies. 



